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On prdsente une approche par dontaines pour le calcul d'gcour 
laments compressibles par des m^thpdes instatiomaaires ou 
pseudo-instationnaires . Cette approche est fondle sur une 
technique gendrale de raccord de deux domaines, dans 
lesquels on rAsout des systSmes hyperboliques qui peuvent 
etre diff brents, I I’aide de relations de compatibilitl 
associies I ces systemes. On montre des exemples d’ appli- 
cation de cette approche au calcul d'dcoulements traasso- 
niques de fluids parfait, et en particulier I I’ajustement 
des chocs. On applique ensuite cette approche au traitement 
d*un probl^me a’ interaction choc-couche limite dans un 
canal transsonique. 


I - INTRODUCTION - 

La mdcanique des fluides numerique s’int^resse S des probllmes de plus en plus 
complexes, en particulier dans le domains adronautique, que ce soit du point de 
vue des configurations geom^triques trait^es ou du point de vue de la moddlisa- 
tion mathdmatique des phdnomenes physiques. Une id€e gin4rale, vers laquelle sem- 
blent converger actuellement de nombreux efforts pour le traitement de problimes 
complexes, et A laquelle d’ailleurs se rattachent des travaux anciens, est cells 
du calcul par domaines. Dans cette approche, le problime du calcul de I’Scoulement 
dans un domains important est dgcompxxs4 en plusieurs sous-prob limes couplds de 
calcul d'icoulements dans des sous-domaines , <S)z y ... dont la reunion est 
Sgale A . Au niveau de chaque sous-probllme , on peut escompter des avantages 
importants de deux ordres : ceux lids i la giometrie et ceux lids au modSle mathd- 
matique ou numerique. 

II sera en effet plus facile de construire, pour chaque sous-domaine, un maillage 
bien adaptd au sous-problime correspondant si les sous -^domaines sont effectivement 
ddconnectds en ce qui conceme leur maillage. Une meilleure adaptation se traduira, 
A prdcision donnee, par qne rdduction du nombre de points de calcul. 

D’autre part, 1 * approche par domaines donne la possibilitd d’utiliser des modllas 
mathdmatiques diffdrents et des mdthodes numdriques diffdrentes dans les divers 
sous-domaines , et done d'utiliser le moddle et la mdthode les mieux adaptds A la 
nature de 1 * dcoulement dans chaque sous-domaine (par exemple, effets visqueiix 
ndgligeables ou non, ou justiciables de 1 ’ approximation de couche limite, dcoule- 
ment irrotationnel, lindarisation I grande distance, etc.). 

Ainsi 1’ approche par domaines permet de tirer profit au mieux des connaissances, 
pricises ou floues, que I’on a a priori sur I'dcoulement I calculer. La difficul- 
td principale que soulive cette approche vient du couplage des sous-prob limes 
entre eux. 



Cette approche peut etre vue egalement sous un autre aspect plus spicifiquemant 
informatique, ^ I’utilisation de calculateurs ^ structure paralllle macrosco- 
pique (multiprocesseurs) , (voir par exemple [1], [2])* Dans ce cas la dicosrposi- 
tion par domaines peut etre recherchee uniquement pour one paralldlisatiou opti- 
male du calcul conduisant ^ une reduction considerable des temps de calcul. 

Cependant le calcul par domaines prS sente des avantages pour le trait ement de 
problSmes complexes m&e sur un calcul at eur convent ionnel (sequential) et seul 
cet aspect est abordS dans le travail pr^sentg ici. 

L* approche par domaines a dtS consider^e dijS depuis quelques anndes pour la r^so*- 
lution de problimes elliptiques (par ex. [3]) ; plus rSceimient elle a it§ mise en 
oeuvre pour des icoulements de fluide visqueux incompressible [4]^ [5], [6], et 
pour des ecoulements transsoniques visqueux [7]. Les travaux de Le Balleur [8] 
sur les mSthodes de coup 1 age fluide p^ar fait- fluide visqueux r assort ant Egalement 
de 1* approche par domaines. Enfin la technique de sous-maillage et de resolution 
par zones [9] qui est mise en oeuvre ici (§ 4.3,7) pour rdaliser un maillage 
suff isamment fin dans la couche visqueuse est un exemple d’ application de cette 
approche au probllme de 1* optimisation du maillage. 

Dans le present article , nous prisentons une technique de calcul par domaines pour 
la determination d’dcoulements compress ib les, visqueux ou non, par des mSthodes 
de type instationnaire explicite. Le raccord des domaines entre eux est effectuS 
d’une fagon syst^matique par la mithode des relations de compatibility, ce qui 
suppose que les effets dissipatifs sont nygligeables au voisinage de la ligne de 
raccord- AprSs I'exposi du principe gyndral de raccord (§ 2), on montre des appli- 
cations en fluide parfait en particulier pour 1 'ajustement des ondes da choc 
(§ 3). 

On prd sente ensuite une methode de riso lution des Equations de Navier-Stokes en 
fluide compressible et sa mise en oeuvre dans une approche par domaines pour le 
trait ement numerique d'un problSme d^ interaction choc-couche limite dans un canal 
transsonique (§ 4) . 

2 - CALCUL PAR DOMAINES POUR DES FROBLEHES HYPERBOLIOUES - 
2.1 - Principe general du traitement des fronti^res - 

On considSre un syst^me d’yqtxations aux diriv^es partielles du premier ordre, 
ecrit sous la forme conservative suivante : 

( 2 - 1 ) ^ Q 

dC 

oD les /V composantes du vecteur U sont les variables indipendantes dScrivant 
completement I’etat du systems physique considgry, t est le temps et Xj des 
coordonnyes d’espace. Les vecteurs ^ , en nombre ygal I la dim^sion de 
I’espace ^considiry , sont des fonctions explicitement connuea des composantes 
du vecteur U . Divers aystSmes telsque (2.1) peuvent etre considyrya en dynamique 
des gaz. Le systime de base est celui des yquations d'Euler qui traduisent les 
lois de conservation de la masse, de la qtaantit€ de mouvement et de I'inergie, 

avec .oaf est la masse volumique, V le vecteur 

vitesse et £ I’ynergie totals spydfique. On peut aussi envisager des systSmes 
simplifiys dyrivys des equations d’ Euler en utilisant des propriytys exactes 
connues de la solution cherchee, ces propriytes poxxvant etre exactes pour des 
solutions instationnaires ou seulement pour des solutions stationnaires. Dans ce 
dernier cas on a affaire S des mythodes de type pseudo- instationnaire dans 
lesquelles le temps apparait comma une variable d*ityration. Le lecteur trouvera 
des exemples et des applications de methodes pseudo-ins tationnaires dans les 
ryfyrences [10], [11], [12] et [15]. 
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A la forme conservative (2.1), nous associons un systime quasi-lin€sire„ Equiva- 
lent pour des solutions continues, que nous Ecrivons da la maniire suivante : 



od et /^' sont des matrices M ^ M dif inies par Aq 
et oil la correspondance entre / et (/ est biunivoque. 


W 



La technique de raccord de domaines que nous prEsentona ici n’est valable que 
pour des systemes hyperboliques dont nous rappelons les propriEtEs essentielles. 
Pour line Etude plus complete des systems hyperboliques, le lecteur pourra consul- 
ter les rEfErences J 13] et [14], Le systeme (2.2) est hyperbolique en temps si, 
pour tout vecteur ^ non nul de , de composantes h , 1* Equation caractE- 
ristique en A : ^ 

(2.3) det ! A fig - ftj j „ 0 

admet^ ^ ]^acines rEelles, distinctes ou non, valeurs propre^ de la matrice 
% ) » et si les t1 vecteurs propres I gauche cC , tels que ! 

(2.4) oC. ( ^ ^ 0 

sont linEairement independants , _constituant ainsi une base de I'espace dans 
lequel on cherche la solution f . 

Dans ces conditions, il existe, pour cheque vecteur propre cC , one combinaison 
linEaire des Equations (2.2) dEfinie par : 


(2.5) 


( K 


df 

dt 



0 


qui ne fait intervenir que 4 dErivEes partielles au lieu de 1 pour le 
systEme initial (2.2). Cette combinaison particuliEre, appelEe relation de compa- 
tibilitE, pent encore s’Ecrire sous la forme : 


( 2 . 6 ) 


ai.f £- i. A fit )i. Jf 


(tj 


oC. L 


oQ les opgrateurs, ... 

reprEsentent des dErivations le long des directions 




vecteurs 


sont contenus dans un plan^ orthogonal au vecteur ^ . II est alors 
clair que si les dErivEes partielles de f suivant les directions sont 

connues E 1* instant Cq , l*Equation (2.6) peut etre considErEe com® un® Equation 
de transport le long de la droite de pente ( V)k ) du plan T ) et permet de 
dE terminer la solution f I 1’ instant tf >^o . Cette interprEtation constitue 

le fondement de la technique gEnErale de traitement des conditions aux limites 
proposEe par Vivi^d et Veuillot dans [11], et dont la technique de raccord mise 
en oeuvre dans le calcul par domines n'est qu^une extension. 


f , 7^^^ j 


et oil les 


</- 1 


On considEre un pc^nt P de la frontilre £ , Eventuellemnt mobile, du domain® 

de calcul ^ ± ^ sst la normale unitaire en P orientEe vers I’eaetErieur du 
domaine W dEsigne la vitesse normale de dEp lac ement du point P : 

W a w . Au point P , le systEme peut etre remplacE par an systEme Equivalent 
composE des reJLations de compatibilitE de la forme (2.6) et dEterminiesen prenant 
pour vecteur ^ la normale extErieure V . Avec ce choix, les dErivEes partiel- 
les le long des directions contenues dans le plan tangent en P E la frontiEre £ 
sont alors connues lorsque les variables ^ sonj: connues dans le domain® . La 
position de W par rapport aux valeurs propres X('p) ordonnEes 4® f aqon 
croissante, par exemple : 

(2.7) 4 .... 4 Ab4 < W 4 Xntif 4 .... 4 Aff 
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OKDntre que, parmi les /V relations de compatibiliti (2.6), m. relatioq.8 
transportent une information de I’extdrieur du domaine de calcul eD vers l*int€“- 
rieur, et ^ ^ m. relations correspondent S un transport dans la direction opposie. 
Les fn premieres relations de compatibility ne peuvent done etre utilisyes et 
doivent etre remp laches par conditions aux limites dont la nature depend 
ggniralement de la physique du problSme considiri. D’oil I'inoncy du principe 
suivant : 

- en un point d’une fronti^re du dotaaine de calcul, les variables sent dyterminies 
I I'aide d*un systSme d’^quationa constituy des /V ^ m relations de compatibility 
associees aux valeurs propres telles que : 

(2.8) A ( v5 j > W 


et de ^ conditions aux limites. 

On notera que, pour le cas limite « U/ , la relation de compatibility corres- 
pondante const itue une Squat ion de transport le long de la frontiSre X et est 
utilisde pour le calcul de la solution. Pour une frontiire fixe ( W * 0 ) et 
dans le cas des systimes pseudo-instationnaires que nous avons dtudiys, le nombre 
m, de conditions aux limites 3. imposer en un point de la frontiSre dypend 
exclusivement de la position du nombre de Mach normal Mp » V. par rapport 

aux valeurs : - I, 0, + 1. En revanche, pour une fronti3re mobile ( ^^(7 ) il 
n'est, en gindral, pas possible de mettre en gvidence j^e jrlgle analogue faisant 
intervenir le nombre de Mach normal relatif Mp a f 1^- v7 , come c*est le 

cas pour les dquations d* Euler exactes. 


La mise en oeuvre pratique du principe pryeddemment gnonci dans une rndthode de 
rgsolution numgrique du systSme (2.1) peut etre rdalisde trSs^ sis^lement sans 
faire intervenir directement les caraetdristiques du plan { ^ , t ) , mais en 
diserdtisant directement la relation de compatibility (2.6) dcrite sous la 
forme : _ 

dans le maillage de base. Aprds diserdtisation, la relation de compatibility 
(2.9) peut etre ecrite de la manidre abrggee suivante 

(7"^' _ (7" 


( 2 . 10 ) 


cC\(. 


dt 


{^r) 


0 


oil y ^ reprdsente une approximation mimerique de • La fagon 

la plus commode de ddfinir cette approximation est d’dcrire le syst^ii (2.1) 
diserdtisd 3 la frontidre en adaptant le schdma aux diffdrences pour tenir compte 
de 1’ absence de points de maillage extdrieurs. Ce systdme diserdtisd se prdsente 
alors sous la forme : 


( 2 . 11 ) 






0 


dt 

prend pour approximation ^ 

..... i . {A . 


dans (2.10), 1* expression 
7^ 


, On notera que les variables (/* n’ont pas d® sens physique 


et on 

diserdtisde 

puisqu’elles n^ vdrifient pas les conditions aux limites. Par dlimination du 
terme dans les expressions (2.10) -et (2,11), on obtient la relation 

trds simple suivante : 

a.i2) 5.\ {u‘" - 0’ ) ^ 0 

ou encore, en employant ies variables / , 


( 2 . 13 ) 


5*. /V 
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Les M^m. relations de compatibility utilis^es en un poin^de la frontiSre sont 
done discrytisges de cette fa<;on en employant la valeur qui risulte de 

I’intygration nnmyrique de (2.1) I la fronti^re. A ces relations on’ adjoint m 
coijjJitions aux limites et l*on risout ce syst^me pour obtenir la valeur de 
I la frontiSre. 


2.2 - Application an raccord de deux sous-"domaines 

On se propose, maintenant, d'appliquer le principe giniral ytabli au paragraphe 
2.1 au cas oil la fronti^re Z est cc^roine 1 deux sous-domaines et . 

Cette frontiSre peut itre une coupure artificiellement iatroduite dans le domaine 
de calcul cD * O , ou bien^encore ime ligne ou surface de discontinuity 
pour les variables fondamentales U ddfinissant, dans le domine ^ , l*€tat 
du systime physique coasidyrd. Si la frontiSre^ est mobile, la vitasse nonaale 
de ddplacement d’un pointy de X est ; ^ ^ 

( a* « f ) . Les noB^res, /7 et /y , des inconnues 

relatives aux sous-d^tsaines et <0^*^ peuvent Itre a priori diff greats. Ces 

inconnues, notdes et , sont lea solutions de deux systimes hyperboli*- 

qxies distincts, de la forme (2.1), auxquels correspondent les systgmes quasi-* 
lindaires suivants : 


(2.14) 


dt ^ dXf 

<£) 

B, £L , Bj ^ ^ 0 

Jm <D 


On notera A et Jl les valeurs propres solutions des deux Equations caracteris- 
tiques, de^la forme (2.^), definies dans les deux sous-domaines 3)^^^ et 
( 9z) ) • Parmi les M (resp, M ) valeurs propres A 

(resp. JL ) ^ ^ (resp. q ) valeurs vdrifient les inggalitgs (2.8) soit : 


4 .... 4 A^^^ < 4 •• 

.{ 




•• « Ji-H 



Dans ces conditions, relations de compatibility doivent gtre e^loyyes dans 
le sous-domaine pour le calcul des M composantes du vecteur / au point P 
de X , et 9 relations sont I utilizer dans le sous -domaine pour le 

calcul des N cotiposantes du vecteur ^ au meme point P . Soit au total 
( ^ ^ ^ ) relations pour dy terminer { ^ ^ M ) inconnues auxqtielles il faut yventuel- 
lement ajouter la vitesse normale de dgplacement du point P ( ou . 

Le p rob lime genyral du raccord des deux sous-domaines et 3)^^^ ne peut Itre 

dyterminy que dans la mesure oCl les ( ^^^ ) relations de compatibility peuvent 
etre compiytyes par ^ relations supp lament aires que nous appe Herons conditions 
de raccord, I'entier i et^t donne par ; 


(2.16) 


avec <f » f lorsque la vitesse VJ est une inconnue et £ » 0 dans le cas 
contraire . 


On peut pryciser le denombrement des relations de compatibility \ utiliser et 
des conditions de raccord dans le cas particulier, taais important, o^ les deux 
systymes (2.14) sont identiques dans les deux sous-domaines et )• 

Nous envisagerons successivement le cas ou X est une coupure artif icielle, puis 
le cas oD X est une ligne ou une surface de discontinuity. 


■:? 

;c 

i 
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-* Cas d*uTie coupure : 


Les conditions de r accord traduisent alors la continuity des M inconnues 
Ui , au point P de la frontiSre L , so it /V relations ; , 

i A 1 ^ .... f M .La continuity de 4/ a pour consequence immediate que les 
valeurs propres ddfinies dans et et les vecteurs propres associes, <i 

et p , sont lies par les relations ; 


(2.17) i 


a) 

b} a Pm. hi 


.... y 


M 


La relation (2.17a) indique que ; CJa . Supposons d’abord que ae 

soit pas ygal I I’une des valeurs propres . La relation (2.17b) montre que 

les p relations de coropatibility S utiliser dans et les (j » rela- 
tions de compatibility I utiliser dans sont liniairement inddpendantes ; 

ces M relations de compatibility inddpendantes permettent done de diterminer 
la solution au point P , la vitesse etant donnye arbitraire. Lorsque la 

vitesse normale de dyplacement du point P est ygale I une valeur propre d*ordre 
r , I cette valeur propre correspondent r relations de compatibility utili- 
sables dans cheque sous-domaine, mais ce sont les memes relations dans lea deux 
sous-domainespuisqu’elles ne font intervenir que des derivies contenues dans le 
plan tangent en P I Z et done identiques par continuity. Au total il reste 
M relations de compatibility indypendantes utilisables. 


-Cas d’une discontinuity : 


Le systyme d’yquations (2.1) dtant sous forme conservative et valable dans^les 
deux sous-domaines et , il est possible de chercher la solution U 

comme solution faible du systSme (2.1) (voir, par exemple, [16]) admettant 
des discontinuitys 2l travers une surface X de position inconnue. Les relations 
de saut associyes au systyme (2.1) s’yerivent : 

(2.18) \>P Fp* ^ IV 

oil 1’ indice supyrieur se rapporte au sous-domaine correspondant i. chacune des deux 
faces de £ . Ces ^ relations scalaires, non lindaires par rapport aux varia- 
bles Ui , constituent les conditions de raccord entre les deux sous-domaines 
et Si la position de (ou de ) par rapport aux 

2 M valeurs propres et est arbitraire, le probiyme du raccord de 

et est, en gyneral, sur ou sous -dyt ermine. En effet, contrairesent A ce 
qui se passe lorsque U est continu A travers £ , la position de par 

rapport aux valeurs propres A| definies dans n’implique, a priori, 

aucime position particuliere de par rapport aux valeurs propres defi- 

nies dans . Dans I’hypothyse oD les relations de saut constituent bien 

AJ relations indypendantes et non dygdnyrees, alors le raccord des deux sous- 
domaines Q,t sera possible si le nombre de relations de compatibility 

utilisyes dans les deux sous-domaines est : /i 1 • 


Cette situation favorable se rencontre lorsque la discontinuity X constitue un 
k-choc , c*est-I-dire lorsqu’il existe un entier ^ tel que les valeurs propres 
A et jc vyrifient la "condition d’entropie" de Lax [17] t 


(2.19) 


< Ai 

AmJ < ^ ^M.hi 
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Dans la sous-domaine on utilise alors } relations de compatibility 

et on en utilise i dans le soua-domaina soit un total de (MiP) , 

Jointes aux Af relations de saut (2-18)* ces (ZMtf ) relations de compatibility 
permettent de calculer les ( inconnues au point P de Z , S savoir : 

les A7 composan^s du vecteur dans la sous-domaine , les Af composan- 

tes du vecteur dans le sous-domaine et la vitesse nonaale de ddplace- 

ment de la discontinuity Z , 

Dans le cas oil Z est une k-discontinuitg de contact , c*est-I~dire lorsqu'il 
existe un entier i tel que le^Valeurs ^ro^res X et vdrifient les 

conditions : 

( 2 . 20 ) 


la conclusion est, en revanche, mo ins yvidenta- En effet, le nombre total de 
relations de compatibility 1 utiliser dans les sous-domaines et est : 

, od r ddsigne I'ordre de la valeur propre A/ (ou )• 

vitesse (ou ) dtant dyterminee par la relation (2.20), le probiyme 

du raccord des deux sous-domaines et «0^^^n'est dont parfaitement dytermind 

que dans la mesure oii, parmi les M relations de saut (2.18),r relations sont 
dygeneryes, c’est-S-dire sont vdrifiyes automatiquement et disparaissent. 

3 ~ APPLICATION AUX ECOULEMENTS DE FLUIDE PARPAIT COMPRESSIBLE - 

La technique gendrale de traitement des frontiSres extemes ou internes ddcrite 
au paragraphe 2 pour des problemes hyperboliques quelconques est dyveloppye ici 
pour les yquations d*Fuler completes et pour quelques systemes pseudo-instation- 
naires . 


3.1 - Systeme des equations d*Euler - 


Ces yquations sont rappelees au paragraphe _4. 2 ,En dcoulement plan, les quatre 
valeurs propres relatives I une direction ^ sont donnyes par les expressions : 

(3.1) Aj ^ d ^ » V. ^ 


'I / ^4 » 

Les relations de conq)atibility s^dcrivent tris simplement sous la forme (2.13) 
en choisissant f ^ i ^ 7 est le vecteur \initaire normal i ^ 

et Vy » ]/ ,rj , et avec ^ [ 0 ^ ^ , 0 , i ) , a 0, 1 ,0) ^ 5^ a 0^0, 1 ) 

^4 M ( 0 , ps , 0 , 1 ) - En ce qui conceme le nond>ra m de condi- 

tions aux limites I imposer sur une fronti^re, mobile ou non, le principe ynoncd 
au paragraphe 2.1 conduit i une discussion simple oil n*intervient que le tiombre 
de Mach normal relatif : 

Afp < - f (frontiSre amont supersonique) : ns 4 

< 0 (fronti^re amont subsonique ou sonique) : m s 3 

0 4t fly < i (frontiSre aval subsonique, ou paroi) : n s 1 

1 ^ ffl^ (frontiere aval supersonique ou sonique) : n a 0 


Considdrons plus en dytail le cas d’une frontiers aval subsonique et le cas d’une 
paroi. Dans ces deux cas, il faut une condition aux limites, et on doit utiliser 
les relations de compatibility assoclies i , Aj , qui s’ycrivent, sous 
la forme (2. 13) : 


.Jifi 


- 








^7 

rp«»;“ I// 


(3.2) 
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et oxlj rappelons-le, I’astSrisque disigne des valetirs fournies direct^ant par 
le schema num^rique S la fronti^re. Pour une frontiSre aval subsoniqua, on uti- 
lise en g^niral ccxoame condition aux limites la donn^e de la pressiom ; connais- 
sant , on tire imediatement et du system (3-2), Pour 

une paroi dont la vitesse de diplacement W est donnie, la condition de glisae- 
ment s*icrit s. W et on determine tr^s facilement 

a I’aide de (3.2). 

La technique de raccord de deux sous-domaines s’ applique au trait^ent d’un choc 
L lorsqu’on risout les Equations d’Fuler completes dans chacun des deux sous- 
domaines, la condition d’entropie de Lax (2.19) Stant satisfaite. Pour fixer les 
idies, on suppose que le fluide traverse Z de vers ; Z est done 

une frontiere aval supersonique de et une frontiere amont subsonique de . 

On utilise les quatre jrelations de compatibility du c6ti de > ce qui signi- 

fie que les T^leurs fournies par le schima sent conservies sans modifica- 
tion, soit m , et on utilise une relation de compatibility (celle 

associee I JL^ , avec les notatiotte du § 2.2) du coti de ; cette relation 

s’icrit : 

(3.3) +(e^j V V 

et forme avec les qmtre relations de Rankine-Hugoniot un systime non liniaire 
aux inconnues et . 

Pour une surface de glissement, la technique de raccord s’ applique ygalement 
sans probl^me car deux relations da saut (en ycoulement plan) sont dygyneryes, 
et il ne subsists que les conditions de continuity de la vitesse normale et de la 
pression. 

Ces techniques sont mises en oeuvre au paragraphs 4, dans un problSme d* interac- 
tion choc-couche limits, pour le raccord de sous-domaines Euler I travers un 
choc et I travers une coupure, ainsi que pour les traitements de la frontiere 
aval et de I’axe de symetrie. 

^•2 - Systlme pseudo-instationnaire ’’H” - 

Ce systyme, interessant pour le calcul des ecoulements stationnaires iso-ynergy- 
tiques, est formy des equations instationnaires exactes de continuity et de 
quantity de mouvement. L’yquation de I’ynergie instationnaire n’est pas utilisde, 
mais elle est remplacye par I’yquation de Bernoulli stationnaire qui dytermine 
explicitement la pression confflse une fonction des variables de base qui sont la 
masse volumique et la quantity de mouvement [10], [11]. 

Ce systyme est hyperbolique, et, en ycoulement plan, les trois valeurs 
propres, qui sont distinctes, s’ecrivent : 

(3.4) A, 3 » + l/-~ ^ ^ f , Ai - Vj 

Avec les variables f coefficients des relations de 

compatibility (2.13) ont pour expression : 

® (f m 

Le nombre de relations de compatibility A utiliser en un point d’une frontiyre 
est : zyro si W < (il faut done imposer toutes les variables) , un si 
Xf 4 W < Xz (il faut deux conditions aux limites) , deux si Xt 4 W < 

(il faut une condition aux limites), et trois si X 3 4 W (pas 
de condition aux limites) . Ces divers cas ne se ramynent pas aux positions du 
nombre de Mach normal relatif, , par rapport A - 1, 0, + 1, lorsque W^O , 

Remarquons que, la methode n’etanc valable que pour le calcul d’ycoulements 
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stationnaires, ^ convergence V! doit tendre vers zSro, et on retrouve une 
classification des^ divers cas en fonction uniquement du nombre de Mach nonaal 
(alggbrique) . 

Les relations de saut S travers un choc correspondent au systeme different 
des relations de Rankine^Hugoniot habitnelles par le fait que c*est I'enthalpie 
to tale (par definition du systhme ”H”) et non l*enthalpie totals relative, qui 
se conserve I travers le choc* La condition de Lax (2*19) s’icrit ; 

(3.5) a) < Aj ! <■ < /^4 

D'aprSs la condition (3.5a), il est possible d*avoir un choc avec un nombre de 
Mach normal relatif amont inf^rieur a I'uniti. L*ltude des relations de choc 
**H” montre qu’I travers un choc v^rifiant (3.5), la pression statique augmente, 
la Vitesse normale diminue, comma c'est le cas avec les relations de Rankine- 
Hugoniot, mais que la masse volumique peut diminuer. L*entropie d’une particule 
fluide traversant le choc peut dgalement diminuer. Bien mtendu ces comport ements 
non physiques disparaissent ^ l*4tat stationnaire. 

Quelques exemples d* applications du syst^me sont prdsentis ci-dessous. II s 
ont Std obtenus par la mdthode explicite de discretisation directs dans le plan 
physique prdsentde au paragraphe 4.3. Le traitement numdrique d'une coupure a 
d’abord etd mis au point pour un dcoulement transsonique, avec choc, dans un 
canal plan. Le maillage utilise est continu I la travers de de la coupure. La 
figure 1 repre sente les courbes iso-nombre de Mach ob tenues I convergence. On 
no ter a que ces courbes ont des tangent es parfaitet^nt continues sur la coupure 
y cotapris' dans la rdgion oQ le choc est capturd. 



Fig. 1 - Profil circulaire dans xxn canal Courbes iso-^ch 

Les autre s essais nximdriques ont portd sur des dcoulements dans une tuydre plane 
utilisde d I’ONERA pour I'dtude expdrimentale de 1 'interaction onde de choc** 
couche limite. Cette gdomdtrie correspond dgalement au calcul prdsentd au paragra- 
phe 4. On montre d'abord un calcul de tuydre compldte, en rdgime amorcd sans choc. 
Le domaine de calcul a dtd sdpard en deux sous-domaines d maillages trds diffe- 
rents (fig. 2a). Le rapport entre la taille des mailles de part et d'autre de la 
coupure est dgal d 4. Etant donnd que I'on ne discrdtise pas d travers la coupure, 
la diffdrence de taille des mailles d travers la coupure peut etre quelconque 
sans qu'il y ait perte de precision thdorique du schdma numdrique. La figure 2b 
reprdsente les courbes iso-nombre de Mach obtenues d convergence. Ce premier 
calcul a permis de ddfinir les grandeurs adrodynamiques dans une section superso- 
nique situde en aval du col sonique. 





b) Courbes iso-Mach ( »»025) 

Fig. 2 - Ecoulement amorc^ sans choc dans la tuySre S8 



U 


q) Courbes iso-^ach ( AM««005) pour p2 ** .67674 pi 
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Pour les deux calculs d’ ^coulements avec choc, le domaine de calcul a iti 
limits au divergent de la tuyere entre la section prdcedemment d^finie *et une 
section situie dans une region od les deux parols de la tuyiire sont pratiquement 
parallSles, ce qui pem^t de justifier I’hypothSse d^une pression statique 
constante dans cette section., Le domaine de calcul c0 a dtS divisi en deux sous- 
domaines c^^^'^et o^^^dont la frontline commune oaobile est le choc Z , ligne de 
maillage particuli^re pour et pour , Le maillage de cheque soua-domaine 
doit done se diformer au cours des iterations pour suivre le diplacement du choc, 
comme dicrit au paragraphs Un premier calcul d'icoulement avec choc a iti 

effectui pour une pression aval p-t » 0,67674 , oil pi dislgne la pression 

d’arret I I’infini anont (fig, 3a). L’ approximation monodixsensionnelle de la loi 
des aires a iti utilisie pour diterminer les conditions initiates, y compris la 
position du choc. On a constati que cette position a peu varii au cours des itira- 
tions. Les risultats de ce premier essai ont servi de conditions initiales I un 
second calcul d'icoulement avec choc, rialisi pour une pression aval plus ilevie: 
fi « 0,68674 pi (fig. 3b). L’onde de choc, et done le maillage, se sont alors 
nettement diplaces au cours des itirations. 

3.3 - Raccords entre systemes pseudo-lnstationnaires diffirents - 

II peut etre intiressant pour le calcul de certains icoulements, d’utiliser une 
mithode met tan t en oeuvre des sys times pseudo-ins tat ionnaires distinct s dans 
deux sous-domaines et * Par example, s’il existe une rigion oU I’icoule- 
ment est irrotationnel, on peut envisager d’utiliser une mithode I enthalpie 
totale et entropie constantes (mithode pseudo-instationnaire '*HSr’ ou ’^2”, 

[10]) dans cette rigion, et de risoudre les iquations d*Euler, ou le systSme 
en dehors. Aucun calcul n'a encore iti effectui dans cette voie, mais il semble 
intiressant de montrer la possibiliti d’un tel raccord, dans le cas, par exemple, 
de la mithode 


La formulation '^S2" correspond au systeme (2.1) avec, en icoulement plan ; 


(3.6) 


U 


<r<i 

<rv 




fM.1 


F, . (<’“f 


oh <r, p et ^ sont des fonctions connues du module V de la vitesse. La fonction 
<r est choisie de telle sorts que le systeme soit hyperbolique. 

Pour une direction ^ donnie, I’iquation caractiristique (2.3) s'icrit : 

(3.7) <tM A" - <r^ ( U ^ ^ \ ^ ^ J ^ 0 

/I - . uSLy 

<TdY 

L’itude de ce trinome en A indique que le systime n’est pas hyperbolique pour 
n’importe quel choix de la fonction (T . En prenant par exemple : 
oh est une vitesse de rifirence donnee, les deux valeurs propres sont rielles 
et ont pour expression : 

( 3 . 8 ) ^ ^1 (f.^j - j , 

Avec les variables / ^ ( V| / Vy / ? les coefficients des relations de compati- 

biliti s’icrivent : 

(3.9) .f- V , eVf(i.^))' 

Traltement d*une coupure **HS2-H* * : 


Considirons par exemple le cas d’un icoulement transsonique avec choc autour 
d'un profil, Le domaine de calcul peut etre divise en deux sous-domaines : un 
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sous-domaine 3)^ od I’^coulement est potential et oil la mdthode **HS2” est utili- 
s€e, et un sous-domaine o^^^^oil l*§coulement peut etre rotacionnel et oS le systS- 
me est risolu. 

Le probleine du raccord ^ travers la coupure Z suppos^e fixe est alors correcte- 
ment posi lorsqne le fluide va de vers et lots qu* on iaipose la continuity 
de I’ycoulement (continuity de la vitesse et de la masse volumique, assurant la 
continuity de I’entropie, soit trois relations). En effet, on dispose de f, rela-“ 
tion de compatibility dans et de f dans ^0^'^ , avec , Cjai si la 

vitesse normale est subsonique, et avec ji^Z ^ qaC si la vitesse nonaala est 
stipersonique. On a done un systime formy de deux relations de compatibility et de 
trois conditions de raccord, pour dd terminer les cinq inconnues (deux dans , 
trois dans 

Traitem^nt d^un choc : 

Soit maintenant un choc mobile Z separant un sous-domaine ^ont I’ycoule- 

ment est irrotationnel et oil I’emploi du systSme *’HS2” est justifiy, et un sous-* 
domaine aval I’ecoulement est rotationnel et oil la mythode est uti- 

lisde. En supposant faible la vitesse de Z , le probiyme est alors bien posy 
en utilisant : 

« les deux relations de compatibility dans , c*est-I-dire en conservant les 
valeurs foumies par le schema ; 

- la relation de compatibility associye I dans <0^^^ ; 

- les trois relations de saut du systime **11”. 

L’hypothSse ** W faible" permet d’assxirer que I'on a bien deux relations de compa- 
tibility utilisables dans . On dispose alors de six relations permettant de 
calculer les deux inconnues dans , les trois inconnues dans <90^^, ainsi 

que W , 

^ “ APPLICATION A UN PROBLEMS D * INTERACTION CHOC-COUCHE LIMITE - 


4*1 - cynyralitya - 

Nous prysentons maintenant une application de I’approche par domainea au traite- 
ment numyrique du probleme de 1 ’interaction choc-couche limite turbulente dans 
un canal en regime transsonique. Ce problime a dejl yty abordy, comme ici, par 
la ry solution des equations de Navier-Stokes moyennes, mais en ne cons idy rant 
qu’un domaine de calcul, le choc ytant capturi dans tout le domaine (par ex. 

[18], [19] et [20]). 

L’approche par domaines permet d’une part de decoupler les probllmes de nsaillage 
dans la rlgion de fluide parfait et dans la region de fluide visqueux, et d' autre 
part d’ajustar le choc sur sa plus grande partie. II en rdsulte une economie de 
points de maillage et une meilleure pr€cision. 

Le canal considery est plan et symytrique. Le domaine de calcul est done limite 
de faqon naturelle par la paroi en has et I’axe de symytrie en haut, et en outre 
de faqon plus empirique par une section amont et une section aval prises assez 
loin de la zone d* interaction. 

Ce domaine (fig. 6) est divisy en trois sous-domaines de la faqon suivante : iine 
coupure horizontale slpare un domaine, vers le bas, qui contient la couche vis- 
queuse et que nous noterons c0vjd’un domaine, vers le haut, oQ les effets visqueux 
sont nygligeafales (domaine 3)p ) , et ce dernier est lui-meme divis€ en deux sous- 
domaines <2}f^ et par I’onde de choc. Autrement dit I’onde de choc est ajustye 
dans <£)p , et captur€e dans , ce qui permet de ne pas resserrer le maillage 
pres du choc dans 0p tout en ayant une precision bien supyrieure i celle que 
donnerait une methode de capture. 
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II serait souhaitable d’ajuster le choc egaleiaant dans ; cette solution est 
€tudi€e raais il faut r^soudre la difficulty qua prysente la foraation progressive 
du choc par coalescence d*ondes de compression. II faut no ter aussi que I'ajuste- 
ment du choc dans 3^ ne dispenser ait pas d* avoir un maillage tr§s resserry dans 
cctte zone afin de bien reprisenter ces ondcs de compression qui rysultent de 
1 ’interaction choc'-couche limite. 

Dans les 2 sous-domaines qui constituent on rysout les Iquations d’ Euler 
completes, et dans le domaine e0y on risout les Equations de Navier-Stokes 
moyennes dycrites au paragraphe 4.2 avec une approximation de couche mince 
(§ 4.3.3). Le choc dans <S)p et la coupure horizontale sent trait^s par la 
technique de raccord de domaines de fluide parfait prdsentye au paragraphe 3.1 
(les effets dissipatifs etant nygligeables au voisinage meme de la coupure du 
cotd de <0^ ) . Le point triple, intersection du choc et de la coupure, ndeessite 
un traitement particulier dlcrit au paragraphe 4.3.5. 

4.2 - Equations et conditions aux limites - 

Nous rappelons tout d’abord la forme habituelle des Equations de Navier-Stokes 
moyennes dans le cas d’une modelisation de type **viscosite tourbillonnaire ** 

(voir par ex. [21] pour une discussion plus complete). 

Les variables de base sont les valeurs moyennes de la masse volumique, p , de 
la density de quantity de mouvement, fv ( V est done la vitesse moyenne pondi- 
rye par la masse), de I’energie interne volumique, pe ( e est done I’inergie 
interne spdcifique moyenne pondyree par la masse) . 

En ycoulement plan et en coordonnees cartesieimes yquations de Navier- 

Stokes moyennes s’ecrivent, sous forme conservative : ^ 

( 4 . 1 ) j. 15. i. AK „ AEl 

_ ()t dK du 

oQ U est le vecteur des variables de base : 

( 4 . 2 ) U ^ CM. , , C ^ ~ 

UL,tr ytant les composantes cartysi^nnes de la vitesse moyenne et £ 
I’energie totale moyenne : £ = e 


Tous les termes dissipatifs sont rassembles au second membre de (4.1) ; les yqua- 
tiona d’Euler pour un fluide parfait s’obtiennent en ygalant le premier i^mbi^ 
y zyro. On a les express ions suivantes pour les termes de fluide parfait 5 


(4.3) 


/> . { , fx£^t/L , cMir , ( ^ / 


oD ^ est la valeur moyenne de la press ion - 
fiques constantes de rapport J , on a ; 

( 4 . 4 ) 

les termes dissipatifs ^ G^. sont donnes 

® ( ^ / / 
s, ( 0 ^ Zxp , , 


Pour un gaz parfait \ chaleurs spyci- 

par les relations : 

^ Zyt ^ ^ Zxv 




i 





( 4 . 5 ) 
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(4.6) 


^‘1 “ ^ 

hi = -T ^ -■§■) 




At ) ck 

f^r ^ {^x. ,d^} 


JUL est la viscositi EDolgculaire, variant avec la tesip^ratura salon la loi de 
Sutherland et JLj est la viscosite "tourbillonnaire" ; est le xmtshxe de 

Prandtl, le noasbre de Prandtl turbulent, pans la phase prSsente de cette 
dtude, consacrie S la mise au point de la aethode de calcul par domaines, on se 
c entente d’un modile algebrique simple sachant qu*un tel modile ne pent etre 
quantitativement correct dans une sone de forte interaction. On utilise le modSle, 
du type longueur de melange, de Michel et al. [22] pour des couches limites non 
dicollees : 

(4.7) JlT = f / 

oil Vf est la composante de vitesse tangentielle ^ la paroi, p la distance nor- 
male t la paroi. La longueur de melange i , rapportde I l^epaisseur / de la 
couche limite, est donnde comme une f one t ion universe lie de p/f : 

(4.8) -L . OASS f/jfJLZ 7 /r. <7,41 

f L(0,0SS S) J 

La fonction , dite fonction correctrice de sous-couche visqueuse, depend du 
frottement total C : C a ( Jl t , par 1 * intermddiaire de la variable 

(4.9) 1 . enp I ^ /(2£ K) j 

La relation (4,7) qui donne JL^ est done implicite puisque F depend de JCr $ 
mais on peut la resoudre numeriquement une fois pour toutes et tabular /A 
comme fonction de la variable Z ^ f /dn j • Dans le cas des Equations 

da Navier-^Stokes , on remplace, dans les expressions ci~dessus, ia 

quantity ^ egaie \ I’oppose du tourbillon. ^ 

En ce qui conceme les conditions aux limites ^ associer aux Equations de Navier- 
Stokes ou aux equations d 'Euler, certaines sont dictSes par la physique ou la 
geometrie du probl^e (conditions I la paroi, condition da S 3 rm^trie), et les 
autres doivent etre choisies de fagon partiellement empirique en s'appuyant sur 
1 'interpretation physique du probl§me et sur la nature mathematique des equations i 
on trouvera une discussion generale de cette question dans [21]. 


A la paroi, on impose I’adh 
nul (paroi adiabatique, 

Sur I’axe de symetrie { V ^ cfe ) les conditions de symetrie sont : tTs 0 . 

Dans la section d’ entree, ou I'ecoulement est supers onique sauf dans un voisinage 
de la paroi, toutes les proprietes de I'ecoulement (c* est-I~dire i7) sont fixees. 


erence du f luide ( U. s tX ^ 
- 0 ). 


0 ) &t \m flux de chaleur 
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Dans la section de sortie oil I’^coulement est subsonique, on impose la valeur 
de la pression conmie on le ferait dans un calcul de fluids parfait* La position 
du choc dans le canal depend essentiellement de cette condition. 


4.3 ^ M^thode numgrique ~ 

4.3.1 ~ Schema de diffgrenceg finies * 

On utilise un schima pridicteur'-correcteur explicits qui constitus uae applica- 
tion du sch^a de MacCormack [23] directement dans le plan physique sn mail lags 
curviligne quslconque* Cs schema a d4j^ €tS prdsenti pour les Equations d’Euler 
[11] et pour les Equations de Navier-Stokes [21 ]^ [24]. Dans le present calcul 
nous avons utilise dans une approximation de type "couche mince ” pour les 
termes dissipatifa qui est introduite directement aur les Equations discrStisdes. 

Rappelons le principe de ce sch&a. Pour touts grandeur <p connue aiix noeuds du 
maillage, on obtient des approximations prdcises au premier ordre des dSriv^es 
^/d£ et ^<P/d^ en un noeud M en supposant que p varie lindairement en 
X 9 y sur un triangle flPQ , oD P et ^ sent deux autres noeuds voisias de M 
(et non al ignis avec M ) ; 

- 4>i^) / ji^p - jiW , « 4 ... 

Ce schima ne fait done appel qu’aux 
coordonnies cartisiennes des noeuds 
du maillage ; la relation qui existe 
entre cette mithode de discretisation 
direct e dans le plan physique et 
une mithode de transformation de 
coordonnies est discutie dans [11]. 

A partir des quatre noeuds voisins 
de soit les points N, S, E, 0 
(fig. 4) , on difinit les opirateurs 
de t^qpe (4.10) suivants : d’une 
part (ir et avec les P^int^. 

M, S, E, et d** autre part i , Sy 
avec les points M, N, 0 (une autre 
variante est ividente) . Le niveau 
de temps itant repiri par I’indica 
%{, SI n. M ) , et la valeur 
pridicteur foumie par le ler 
pas itant notie %fi , le schima 

t/'~' - J*- it[S^{F,.6,]'fSf{F^- Sj] 

u"’ , ±^u\ u^’-at[S,(F-?,r\j(F,.^r]j 

Pour discritiser les dii^vies premiires de qui apparaissent dans les 

termes dissipatifs , on utilisera les opirateurs opposis I ceux utilisis 


a’icrit ; 


(4.11) 


i>J 



Fig. 4 - Discritisation directs 
dans le plan physique 
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pour , c*est~S*“dire Jr , au ler pas et ^ fSu au 

second pas. Le schema est du second ordre en espdc e e t en tei^s si le maillage 
est suffisamment r€gulier (c*est-I-dire si ] t MO j s. Of^^) et 
n o(m^) ) . 

4.3.2 - ■“ 

Deux i tapes supplSmentaires de calcul sont ^ventuelle®ant effectuies I la suite du 
correcteur (4.11.b) pour introduire une viscositd artificielle. La pr^iire dtape, 
ndcessaire I la capture des chocs, utilise la viscositd artificielle non lindaire 
de Lerat et Sid§s [25]/ mais restreinte S une seuie direction, celle des lignes 
de maillage du type OM£ (^ack) qui suivent la direction principale de 
1’ Seoul ement : 

U) u”"'" f 

oO JX est de 1 ’ordre de 1’ inverse du nombre total de mailles dans la direction 
des lignes J. s ck , et oil <£ est pris de I’ordre de I’unitd. 

La deuxi^me dtape introduit une viscosity artificielle liniaire du 4e ordre [26], 
permettant d’ amort ir les oscillations de courte longueur d’onde sans perturber 
lea zones ^ fort gradient de la solution : 


(4.13) 


1 ^££' ^OO' 

- -3 

- 3 

j 




— 

' 


jittiji) 


- ] 



- 5 (Jfifs 

i 


oD C/ , /V ^ ^ sont les points de maillage respectivement voisins de 

E , 0 , N f S indiquds sur la figure 4^ et oil les coefficients , Ei sont 
pris de I’ordre de 


4.3.3 - App roximation de co u che mince - 


II est possible de simplifier les calculs en introduisant une approximation de 
type "couche mince" direct ement dans les Equations diserdtisdes. Les lignes de 
maillage du type OME dtant orientdes salon la direction principale de l’dcoule«* 
meat, I’hypothdse de coxache mince revieat I considdrer que les termes dissipatifa 
au second membre de (4.1) sont essentiellement determine par les gradients 
transversaux (c ’est-d-dire dans la direction de StlN ) de £ . Si, dans 

les formules cdndrales (4.10) on convient que P reprdsente soit M soit S , 
et Q soit £ soit 0 , ces fonmiles pourront etre simplifides en icrivant : 


(4.14) 


i^x<p)ls^ = - / I 

^)m ^ ^MP / ^ 


et ce, uniquement pour les termes dissipatifs. Les Equations d’Euler resteut 
done entidr ement contenues dans les equations diserdtisdes, et celles*ci restent 
dga lament valables dans , que ce soit dans la couche visqueuse ou I 

I’extdrieur. 


4.3.4 - Calcul en m ail lage varia bl e - 

Le schdma a dtd ddfini ci-dessus en supposant le maillage fixe. Or le calcul est 
effectud dans un maillage variable car les lignes transversales se ddplacent 
pour suivre le mouvement du choc, aussi bien dans que dans «0y - 

Repdrons les noeuds du maillage de faqon habituelle par un couple d* indices ) 
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s4par€ment pour et pour , Las lignes J&cfe , obtenues par interpola- 
tion en ^ , sont fixes, Dans 3)p le choc esc ajuste ; c*est une li^e de 

maillage mobile, d* indice i^ic fix^, oui coune l* interface entre et 
au Doint triple A (fig. 5) . 


CHQQ 



Fig. 5 - Maillage au voisinage du point triple. 


Dans on prend les lignes <*ck verticales (les formulas (4*10) sont 
simplifides en consequence), et la verticale issue de ^ , x ^ Xp (t) , est 

choisie consne une lizne de maillaae . d’ indice fixS ; ainsi le 

resserrement des lisnes isck qu*on introduit autour de pour ime bonne 

representation de la zone d* interaction, suit automatiquement le d^placemant 
du choc. 


Pour appliquer le schema en maillagg variable on pourrait^ comiae dans [24], 
remplacer dans (4.1) par ^^/dz - Kr reprS sente 

la dSrivde en temps en un point de maillage donnS de vitesse 57 (de composantes 
Kt / ^ ) et intdgrer p^r rapport A Z . Nous utilisons ici une autre techniaue 
consistant I calculer sur un maillaae fixe d^fini ll 1* instant . comme 

dScrit dans les paragraphes precedents, puis A projeter cette solution sur 

un noweau maillage defini A 1’ instant 


Soit M , E et 0 les positions A I’instant ^s points de maillage 

respectivement. La vitesse k4‘/ ^ d^placemant du point 

(^J) est determinee par les 2 conditions suivantes : a) Wc,j est parallSle ^ 0£ ^ 
b) la composante est obtenue par interpolation lin^aire en X entre le 

choc dans oDp , ou la ligne £ a dans , at la frontiire fixe amont ou aval. 
On calcule alors la position du point (^,/) 1 I’instant , ainsi 

que la valeur en M' de toute grandeur <p connue dans le maillage de 1 ’instant 
, par un d^veloppement de Taylor au ler ordre : 


(4.15) 



K-/ = X at 0£ 

6(M) f Wm . <p 

i X ilt [ (p(£) . j>(0)J 
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En pratique la vitesse de dgplacement du maillage reste trSs faible et le 
dSplacement M'M reste toujours tr^s inf€rieur ^ la dimension de la maille M£ 
ou W . 

4.3.5 - Tr ait ement du point triple - 

Le point triple ^ nicessite un traitement particulier assurant le passage de 
la m^thode d'ajusteim^nt I la methode de capture pour le choc, et tenant compte 
de c@ que les lignes i a eft sont beaucoup plus resserr€es dans que dans 
au voisinage de ce point* 

Le schema de la figure 5 montre le detail du maillage au voisinage du point 
triple (les interfaces entre sous-domaines sont didoubldes pour bien distinguer 
les points de maillage) . 

Le_sch&a numdrique n'est pas appliqui aux points /) et > et les valeurs de 
en ces points, qui interviennent dans les relations de cos^atibilitS (2.12) 
sont obtenues par extrapolation parabolique le long de 1 'interface a ch) 

S partir des points amont pour ^ , et I partir des points aval pour .La 
vitesse du choc qui r^sulte du traitement des noints du choc n'est nas retenua 
en 4 pour calculer le diplacement de ce point ; en effet afin de relier le choc 
ajusti dans au choc capture dans , on dit ermine le point choc comme 
itant i 1' intersection de 1 'interface avec le segment BS, S itant le point Soni- 
que sur la li^e J » ck insaidiatement an-dessous de la coupure. Pour les points 

de maillage situis^ sur la coupure du coti de et entre ^ et £ , on deter- 
mine les inconnues U par interpolation liniaire entre et d’une part, 

et entre et D d’ autre part, en prenant U^/r { U/i ) / Z ; en £ et F 

ainsi qu'en C et £? , on applique le traitement de la coupure dicrit au paragra- 
phe 2.2.L’italen^nt du choc k la coupure est ainsi limiti I 2 mailles de • 

4.3.6 - Tra itement d es frontieres - 

a) Paroi s le schima (4.11) n'est pas appliqui ^ la paroi f ) , imsis seulement 

'k partir de la ligne au-dessus (j,aZ) . L'icoulement I la paroi est obtenu k 
partir des conditions aux limites ( mx a 0 , ^ 0 ) , et i I'aide de 

1' Equation de quantiti de isouvement normale qui foumit en fonetion des 

termes dissipatifs. Pour des icoulements 1 grand nombre de Reynolds, et pour une 
paroi k faible c our bur e, on peut prendre « 0 . 

Compte tenu du choix des opirateurs ^ qui interviennent dans (4.11), 

le calcul du correct eur sur Z ne fait appel k aes valeurs pridicteurs sur la 
paroi que pour les termes dissipatifs, et il n’est pas ndeessaire de dSfinir de 
valeurs prSdicteurs pour ^ et f * 

Les formules aux differences, d^centrSes en p , utilis€es pour discr^tiser 

et ^ is paroi, conduisent k das syst^es tridiagonaux que I'on 

rSsout ais^ment par la methode de double balayage. 

b) Front i^re amont ; toutes les grandeurs sont donn^^s sur la frontiers amont 
(ia f j , On y calcule neanmoins \m prSdicteur (J^*^ qui est n^cessaire au 

calcul du correcteur sur la colonne suivante ( ixZ) 

c) Axe de symetrie : e'est une frontiSre du domaine . En fluide parfait, il 
est Equivalent de traiter cette frontiire comme une paroi plane, ce que I'on fait 
par la mEthode decrite au paragraphs 3.1. 

d) FrontiEre aval ; la pression est imposEe, et on utilise la mEthode basEe sur 
les relations de compatibilitE decrite au paragraphe 3.1., aussi bien dans 3p 
que dans Sy , Nature! lament dans , les valeurs sont cel les foumies 
par le schEma pour les equations de Navier-Stokes , et non pour les Equations 

d ’ Eul er . 
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' 4, 3, 7 - 

Un maillage extremement fin est n€cessaire pris de la paroi pour bieti repr^sen- 
ter les variations trSs rapides des proprietes moyennes de I'^coulement dans la 
couche visqueuse turbulente, 

Ce resserrement de maillage est rialis^ en comb inant deux techniques ; a) une 
contraction r^guliSre du maillage de base d^finie par une loi exponentielle : 

! (p' ;/o(^)) = (e^^-1 ) /(e^-1 ) 

oh paroi et ^ coupure, le maillage de base itant, 

par aSf tuition, h pas constant dy s et b) une technique de sous-maillage par 
zones [9]^ [24], dans laquelle la maille en / est divisSe par deux d’une zone 
2L la suivante en allant vers la paroi. Cette technique de dichotomie permet un 
raffinement rapide du maillage, tout en conservant un maillage rdgulier dans 
cheque zone. Le raccord des zones est as surd par recouvrement, et un systSme 
spdcial de balayage en temps des zones est mis en oeuvre de faqon I utiliser 
dans cheque zone un pas de temps maximum propre I cette zone. 

La figure 7 montre le sous-^illage realisd avec cinq zones pour 1* application 
mimdrique prdsentie au paragraphe 4.4. 

4.4 - Rdsultats - 

Les conditions du probleme d ’ interaction choc-couche limite traitd ici sont 
celles d*une experience rdalisde dans la soufflerie SB de I'OHEIA et ddcrite dans 
[27]. Les conditions gendratrices dtaient les suivantes : pression gdndratrice 

» 0,96 bar ; tempdrature gdndratrice « 300®K.Le nombre de Reynolds calcu- 
Id h partir de la hauteur d de la section aval ( d « 50 m, demi*-hauteur du 
canal), de la vitesse du son gdndratrice et de la viscositd cindmatique dans 

les conditions gdndratrices, est dgal S : « 1,039.10^. Le noM>re de Prandtl 

turbulent vaut 0,9 ; le nombre de Prandtl est dgal I 0,725 et le rapport des 
chaleurs spdcifiques k 1,4. Enfin, le nombre de Mach h la frontidre de la couche 
limite et h 1* origin© de 1* interaction vaut environ 1,3. 

Le domain© de calcul, qui s'dtend de » 113 mm I X » 238 mm, et les maillages 
(h convergence)^^^sont reprdsentds sur la figure 6. Le maillage comporte (22 x 13) 
points dans t£)p et (30 x 13) points dans cDp^ ; le maillage de base d© <s2y com- 
prend (61 x 21) points. 

Les conditions sur la frontidre amont Zf ont dt§ ddtermindes 1 partir d*xm calcul 
prdalable, par la meme md thode, effectue dans la partie amont du canal reprdsen- 
tde sur la figure 2. a, pour un rdgime amorcd sans choc. Les conditions initiales 
sont ddfinies dans J)p en orenant une solution monodimens iotmelle corresnondant d 
une loi de section du canal tenant compte d’une dpaisseur de ddplacemant ddduite 
de I’expdrience ; la pression sur la frontiSre aval est ddterminde de telle 
sorte que le choc de la solution monodimens ionne lie se place I 1 'abscissa du 
choc expdrimentale. Les conditions initiales dans sont ddduites des conditions 
amont (cx^l )et des conditions sur la coupure ( ) par les fortmiles : 

(^J ) /^ ^ / 4* pour <f> u u , ft, , H . Sur la frontidre 

aval, la pression est maintenue const ante, dgale I sa valeur initial e dans 
( ^2 ® 0,6625 f,i ). Cette valeur diffdre de la valeur expdrimentale ( « 0,6357^^4^ ) 

par suite des effets tridimensionnels dus d la prdsetice, dans 1' experience, de ‘ 
couches -limites lateralcs. 

Dans le moddle de turbulence, on utilise une loi d* dpaisseur S expdrimentale, 
calee sur le minimum de pression d la paroi. 

Les risultats prdsentds correspondent d un noirfjre d' iterations dgal d 3000, soit 
un temps physique ecoule egal d 3,59 4 /^i « 5 , 16. lO'^^sec. En ddbut de calcul, le 




Fig. 6 ** Oomaines de calcul et maillages. 



Fig *7 Sous-<^aiirage par zones. 



a) Interferomdtrie 


b) Courbes isochores ( 




.02 a ) 



c) Courbes iso-Mach ( AM = .02) d) Courbes isobares ( A|u «■ .01 -|v ) 

Fig* 8 - Ecoulement dans la region d' interaction choc-couche limite 





21 



f20 f4Q Md f0Q 2M xfmm) 

U) ft) fc) (d) 


Fig. 9 - Distributions calcul^es de la preaaion pari€tala, de 



Fig. 10 - Profils de 



22 


choc se diplace vers l*aval, puis revient vers l*aaK)nt et repasse par sa position 
expirimentale (au niveau de I’axe de sym^trie) S 3000 iterations. La solution 
presentie ne correspond done pas I une position d'equilibre du choc, mais on pent 
la considlrer cosane etant en evolution quasi-stationnaire, car la vitesse du 
choc est faible (environ 0,023 ). 

Le module de turbulence utilise ne pouvant etre quantitatives^at correct dans la 
zone d’ interaction, on ne prdsente pas ici de comparaison ddtaillle avec les 
rdsultats expdrin^ntaux. La comparaison est cependant tout I fait satisfaisante 
qualitativement, comma on pent le voir sur les figures 8a et 8b, qui montrent 
des courbes isochores, respectivement expdriaentales (obtenues par interfdromd- 
trie holographiaxie) et calcul€es. Les fibres 8c et 8d representent les courbes 
iso-Tiombre de Mach et las courbes isobares calculdes. On notera sur la figure 8c 
I’Soaississes^nt rapide de la couche visqueuse, provoque par 1 ’ interaction avec 
le choc (une partie seuleaent du domaine de calcul est reprisentee sur ces 
figures) . 

Les variations de f/^ , od % est la pression pari€tale, de l^dpaisseur de 

dgplacement 5^ et du param^tre de forme Hi (rapport de I’Spaisseur de ddplace- 
ment "incompressible’* a I’^paisseur de quantity de mouvement "incompressible”) 
en fonction de X sont reprdsentees sur la figure 9. La solution obtenue comporte 
un Idger ddcollement, d'^paisseur faible, au plus ^gale I 0,15 mm, sur vme 
distance d’ environ 25 nam. Les points de ddcollement D et de recollement R sont 
indiquis sur la figure 9. Dans Inexperience, 1* exploration de la couche limite 
par vdlocimetrie laser a montre que celle-ci est tr^s proche du decollement. 

Des profils de vitesse longitudinals u, sont reprgsentes sur la figure 10, a deux 
Schelles differentes, et pour les abscisses indiquies sur la figure 9 : un profil 
(a) a I’amont du choc, un profil (b) dans le choc, un profil (c) dans le dicolle- 
ment et un profil (d) recolle. Les points de maillage sont port^s sur ces profils. 

Les calculs ont etS effectu€s sur un CYBER 170*750. et le temps de calcul pour 
les 3000 iterations est d’ environ trois heures. Neanmoins. ce temna ne doit etre 
pris cue comma un ordre de ^scrandeur ; il est tr^s sup^rieur au temps de calcul 
qui pourrait etre obtenu avec une optimisation du code numdrique. En effet, le 
calcul presents a dti mene rapidement en associant deux prograiomes existants, 

I’xm pour , 1* autre pour , progratsmes qui avaient Std r^alisis indSpen* 
damment l*un de 1’ autre. 


Cette Stude a et€ effectude avec le soutien financier de la Direction des Recher* 
ches Etudes et Techniques, et du Service Technique des Programmes ASronautiques . 
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